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Proc 

Description: 

The invention is concerned with a procedure to assess deformation during which they are 
determined. Subject of the invention is a precise examination how an influx of energy due to loading 
causes a deformation. 



Furthermore, the invention concerns a mechanism which is suited for the implementation of 
the procedure. 



Several procedures are known by which the deformation of a body due to action of forces is 
assessed. 



The known procedures are based on the method to represent the body through a grid of 
nodes; for every grid node the change of the relative position due to the action of forces is assessed. 



These procedures are rather cumbersome for complex body shapes and especially for very 
precise solutions at small relative scale. For example, the modelling of deformation of a vehicle due 
to a crash requires an extremely high calculation power. 



Although the known procedures employ mathematical models that are in principle fairly 
simple - the employed mathematical models are essentially based on the work of Euler in the 18*^ 
century - there are no procedures known so far that differ fundamentally from the procedures based 
on Euler. 

The inventor has developed a theoretical approach to the theory of deformation which are 
presented in a number of freely available articles. 



The first of these articles, "An approach to deformation theory based on Boyle's law. I. The 
concept, and kinematics of pure shear" offers a concept to assess the kinematics of deformation of a 



2 



solid based on an adaptation of the ideal equation of state PV = nRT. This theory makes it possible to 
consider defonmatlon involving both isotropic and anisotropic materials and/or boundary conditions. 

Further developments to the kinematics of deformation have been presented in the article "An 
5 approach to deformation theory based on Boyle's law. II. Kinematics of simple shear, and some 
energetic considerations". 

The necessity to consider a full deformation in three dimensions even for very simple 
geometric setups is explained in the article "An approach to deformation theory based on Boyle's law, 
10 III. Three-dimensional properties of plane-strain deformation: origin of conjugate joint sets, sheath 
folds in plastic shear zones, and turbulence in viscous flow". 

Subject of the invention is the task to develop a general method to the point where it delivers 
solutions to specific tasks as quickly and reliable as possible. 

15 

In the sense of this invention the problem is solved by the assumption that a point Q within a 
body of solid is thought to serve as a reference point which is separated in all directions from 
corresponding points P through a unit distance. The points P serve as direction indicators relative to 
Q; together they form a surface which has spherical shape under ideal conditions. The change of 

20 location of all P relative to Q as a function of the acting forces describe the deformation at the point Q. 
In the traditional procedure a grid of nodes is pre-defined, and all points serve simultaneously as 
reference points and as correspondence points. Therefore the accuracy of the procedure is 
dependent on the number of point pairs included in the consideration which may be very high, but 
which is necessarily finite, and thus there are gaps. In contrast, in the procedure proposed here any 

25 point can serve as reference point Qi, but the surface of all respective points Pi is then defined; if the 
physical relations at Qi are considered, it is irrelevant that a point Pi may serve itself as a point Q2. 
This way there is an infinite number of P for every Q, and within the defined domain of validity (i.e. the 
shape of the body of solid) there is an infinite number of Q, which makes it therefore possible to apply 
the method as a homogeneous differential approach. 



Previous method 


Method presented here 


The grid of nodes needs to be pre-defined 
(Fig.01). 


Any point may serve as point Q which is 
therefore by definition surrounded by a surface of 
all points P (Fig.01). 


The number of nodes may be arbitrarily large, but 
it is finite. 


The number of points P for every Q is infinite. 


The number of connecting lines between pairs of 
points grows with the number of points 
considered, but it is discrete, and every pair 
requires its own mathematical modelling. 


The number of connecting lines between one Q 
and all corresponding P is infinite. Their 
mathematical treatment therefore permits a 
continuous field method. 


The distance between any two points is finite. 


The distance Q->P is scale-independent. 


Two arbitrary nodes in the grid may serve as 
starting point and end point of a given distance. 


The distance fo is the distance between a Q and a 
corresponding P, but not between two Q. There is 
no relation between the surface point Pi of a norm 
volume centered at Qi and another Q2, (not even 
if Pi is simultaneously some other P2,) except if a 
Pi happens to coincide with a Q2. ' 


Deformation is modelled through a change of this 
distance in length and/or direction relative to one 
node or the other, and relative to an external 
reference frame. 


Deformation is modelled by means of r relative, to 
its own reference point Q only, and only indirectly 
relative to an external reference frame. . 


The modelling of deformation is constrained to the 
X discrete connecting lines of all pairs of points in 
the grid. 


The modelling of deformation is homogeneous- 
continuous over all P relative to their respective 
Q, and homogenous-continuous over all Q 
relative to an external reference frame. 


Deformation is understood as a mathematical 
mapping of mass points from the undeformed into 
the deformed configuration. To prevent the 
mapping of two mass points upon one another it 
is necessary to apply the compatibility equations 
apparatus. 


Deformation is understood as an - energetic 
change of state. Two different points P and Q 
cannot coincide because the work equation is 
logarithmic, thus infinite work would be required to 
let two points coincide. The mathematical 
apparatus of the compatibility equations is 
therefore not necessary. 



Symbolic Terminology 



T, M, F tensors f, m, r, s, v vectors 

I identity tensor n.t unit vectors 

Ty, F;7 tensor components n vector components 

P, V, H, U scalars. state functions f. r vector magnitudes 

P, Q points a, y, 0 angles 



The forces which displace the surface of all P consist of the sum of external forces text, the 
forces exerted by the system upon the surrounding m^^su and the surface-bonding forces mA which do 
not do work themselves (they are constraining forces), but they make work done by stretching 
components of fext and the shear components ms(syst) and fs(ext) possible: 



fext + >Tlsyst + mA - 0. 



(Eqn.01) 
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The new elements in this method to consider loading and defomiation are 

the introduction of a region in space of fixed extent (the system) which is separated from its 
environment, and with which it interacts mechanically, and whose physical properties are scale- 
5 independent (independent of the size of the system), 

the explicit separation of forces fext which are controlled and exerted by the surrounding upon 
the system, and the forces mgyst which are exerted by the system upon the surrounding, 

the inclusion of the existence of bonding forces mA (permanent within a solid, transient in a 
fluid) which considerably constrain the freedom of the system to react upon external forces, and 
10 the assessment of the energetic effect of shear forces (the work done by shear forces) and 

the transformation of the shear components fs(ext) and ms(syst) of the vector fields fext and msyst into 
normal pomponents parallel to the main extending direction. 

The invention thus consists of a method to assess deformation of a body whereby it is the 
main concem how the action of forces affect the body. It is done in such a way that the equilibrium 
15 between internal and external forces is established as a function of material properties, external 

boundary conditions, and the existence of bonding forces between system and surrounding, both for a 
representative part of the body, and for the body as a whole.'The material properties must be known ' 
as a function of the location Q(X|). Likewise, the configuration of the external boundary conditions 
(properties of the external force field) must be known; 

20 

The entire work done by the external vector field is determined through the partitioning of 
normal forces and shear forces according to the equation 

r * f = ^|rxf|' +|r.f|' = |r||f | = const (Eqn.1 8) 

25 

The radius r is a mechanical lever. Within a time period that is short relative to the material- 
dependent diffusivity or the relaxation time, the method is also applicable to fluids (liquids and gases). 
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The total displacement is determined through a work equation that must be derived from the 
equation of state applicable to the material under consideration. In the simplest case the work 
equation is 



5 

jfdr = c|— = clnr . (Eqn.20) 
which is derived from the ideal gas equation of state PV = nRT. 

5 In the following description the invention is explained using two-dimensional examples. An 

expansion to three dimensions can be achieved without difficulty. 

Fig.01 Upper panel: Known method (finite element method): a grid of nodes is pre-defined. 
The deformation at the point Q is determined by considering change of length and orientation of all 
10 distances whose one end is Q. The respective outside points are the corresponding points P. The 
finite element method is named such because the distance between any two nodes Is finite. Lower 
panel: Any point can serve as reference point Q. and all corresponding points P form a closed surface 
around Q. All relevant variables - material properties, boundary conditions, deformation - can be 
presented as differentiable functions of location /(Q). 

15 

Fig.02 Compressibility data for the alkali elements. External load expressed in multiples of 
the internal pressure. Data recalculated after Bridgman (1958, p.180). 

Fig.03 Relation of fop to fdev- Vertical: vector magnitude, horizontal: angular distance. The left 
20 and right side represent the contracting and extending eigendirection (c and e) of an anisotropic field. 
The eigendirections are orthogonal and simultaneously parallel to the principal axes of an ellipsoid 
only for an orthogonal field. Fig.03 shows relation of minimal and maximal external forces relative to 
the operative field. Gravitational forces are not considered. Confining pressure and operative field are 
isotropic-hydrostatic. 

25 

Fig.04 (a) Deviatoric force field for pure shear. fi and f2 are perpendicular to one another 
and parallel to the coordinates and the eigendirections of the field. The same figure represents the 
displacement field or flow field if the arrows stand for displacements, (b) Orthorhombtc entire force 
field, schematically, mf, material force, fop: hydrostatic operative pressure, f,*: deviatoric force field. 
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Fig.05 Force magnitudes along the surface of a unit body for pure shear. Continuous line: fn, 
dashed line: fs. Their sum is unity at all points. 

Fig.06 Dilatancy effect of a shear force. Circular line: surface of a thermodynamic system. 
5- Tangential force vector (dashed) will displace material point P from P to P*; the distance from the 
origin increases. Point coordinate Xi will be stretched by horizontal component of f to Xi\ X2 will be 
shortened by vertical component of f to X2'. Xf. external coordinates parallel to the eigendirections. 

Fig.07 Outer bold line: surface of a unit body before loading. Inner thin line: unit body after 
10 hydrostatic loading due to an operative field, resulting in volume decrease. Inner elliptic line: change 
of shape of the loaded body only due to normal components of the deviatoric field; the volume effect 
is zero. Outer elliptic line: finite deformation after consideration of the work done by shear forces. Area 
between the two elliptic lines: volume increase due to the dilatancy effect. Xf. external coordinates: 
parallel to the eigendirections. 

15 

Fig.08 External force field for simple shear with unit body, (a) The force field consists only of 
forces parallel to Xi. (b) The unit body is fixed in space and cannot rotate due to bonds between body 
(system) and surrounding. P: point of action of a force vector, O: coordinate origin, 9: angle between 
position vector of P and internal coordinate Xi. 

20 

Fig.09 Forces on the surface of a volume element in a continuum subjected to a force field 
as shown in Fig.08. Average force fav acting on its point of action Pk on the surface of the system 
during external conditions of simple shear, fav with magnitude sin k decomposes into normal 
component fn and shear component fg. The latter is balanced by surface bonding forces. fs(K) is 
25 subtracted from all shear forces at all points, resulting in sinistral shear forces between P^ and 0. X is 
the angle of a coordinate transfonmation which results from the imbalance between the shear forces 
over the surface of the system. The equilibrium condition f f x r d9 = 0 is observed. 
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Fig. 10 Thermodynamic system, resp. volume element with unit mass and elliptic shape 
which is in equilibrium with a monoclinic extemal force field. Acute angle between eigendirections c 
and e = 2k. 

5 Fig.1 1 Vector magnitudes in simple shear, (a) in untransformed space, (b) in transformed 

space. Continuous line: normal component fn(dev) (eqn.44), long-dashed line: | f<jev x I (eqn.42), dot- 
dashed line: fs (eqn.45, see text). ±R. ±P : orientations perpendicular to R- and P-surface, cf. Fig. 12. 

Fig.12: Geometry of simple shear, (a) Calculated geometric properties: X/. Xr, coordinates, c: 
10 contracting eigendirection, e, extending eigendirection. ±R: direction from the origin to the point R and 
perpendicular to the R(iedel) plane; ±P: direction from the origin to the P-point, and perpendicular to 
the P-plane. (R, P: terminology cf. Berthe et al. 1979.) (b) observed geometric properties (Weber & 
Juckenack 1990). The only major difference between (a) and (b) is the orientation of the grain shape 
foliation. It is assumed, however, that the extending eigendirection should be correlated with the 
1 5 crystal lattice fabric which is commonly a little lower. 

Fig.1 3 (a) Displacement field for simple shear. Although the external force field does not 
have a component perpendicular to X2. such a component does exist in the effective field f^ev and in 
the displacement field, (b) Natural SC-fabric: mylonite in lower greenschist facies, Insubric Line, Sesia 

20 Zone, Val Strona, Italy. Long dimension of picture is parallel to bulk foliation which serves as external 
reference frame. The shear planes inclined gently to the right are the C-planes. The layering dipping 
shallowly to the left is the S-plane. Considerable sinistral whole-body rotation between C- 
discontinuities which themselves have dextral sense of shear. Bulk sense of shear is dextral. (c) 
Viscous simple shear in subrecent obsidian flow, Lipari Island, Italy. Upper layer consisting of black 

25 glass with bubbles was softer and shows drag features. Lower layer cohsists of partly crystallized 
material and reacted brittle. Drag in upper layer and joint orientation in lower layer both indicate 
dextral bulk shear. 

Fig. 14 Top: rectangular sample of solid between two anvils with the location of points for 
30 which the boundary conditions are shown below. Center: boundary conditions as a function of location 
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within the body. Circles stand for (from inside out): operative field (isotropic); normal component of 
deviatoric field; shear dilatancy contributed by the system; volume effect due to extemal boundary 
conditions. The latter is directed inward in (1) due to the friction effect along the sample-anvil contact, 
in (2) It is directed outward due to the shear dilatancy effect contributed by the surrounding, and zero 
5 in (3) for points on the boundary between solid and freespace. Bottom: presentation of the expected 
effect; white circles: unloaded state, black ellipses: loaded state. 



lower panel are mathematical artefacts due to summation over non-infinite coefficients n. 

10 

.Fig.16 Dimensions 10x5 unit distances. Wavy structure along the horizontal boundaries in 
lower panel are mathematical artefacts due to summation over non-infinite coefficients n. 

Fjg.1 7 Dimensions 10x10 unit distances. 



Fig.15 



Dimensions 10x2 unit distances. Wavy structure along the horizontal boundaries in 
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Fig.18 



Dimensions 10x30 unit distances. 



Model calculations are here restricted to two-dimensional cases for simplicity. In this 



. application equations are commonly given in vector notation and in algebraic form, e.g. 
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J f X r de = J sine cosG d0 



The invention comprises four parts: 



25 



the definition of a scale-independent unit volume with a shape that fulfils the equilibrium 



conditions, 



the definition of two force fields which are in equilibrium with one another, 



the determination of a material-independent equation of state for solids, 



the transformation of the equation of state from scalar form into vector form. 
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Definition of the system 
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The purpose of the system is the definition of a reference body in the unloaded state which 
represents the material properties and other mechanical properties. Due to loading, work is done 
upon this system which is scale-independent in the sense that similar work is done upon a small or a 
5 large system relative to the mass involved for identical respective deformations. Thus it-is ascertained 
that it is possible to calculate with a finite radius or standard volume although their actual numerical 
value is irrelevant. 

Any mass is associated with a potential. In the case considered here the internal energy of a 
10 thermodynamic system is considered as its potential. In a given energetic state it is proportional to its 
mass, dU oc dn, where n is the number of mol (the thermodynamic mass is dimensionless). The 
energy density per unit mass is then dU/dV = P, in integrated form it is UN = P, such that the mass is 
proportional to both U and V. UN = P is then the energy density, or the thermodynamic potential per 
unit mass. The latter is then finite. 

15 

In the reference state Uq we have PV = nRT. An infinitesimal change of state requires the 
additional energy dU, Of course, the nature of this dU at constant n is very different from that of dU oc 
dn above. In order to consider infinitesimal changes of state, a system of finite mass and volume is a 
prerequisite - a thermodynamic system - which is associated with its self-potential. 

20 

In order to define the work done on a point Po, a reference point Q must be chosen such that 
Q ^ Pg. The choice is arbitrary and may be defined to the convenience of the problem (Kellogg 1929, 
p.53, 63; Boas 1983). The distance Q->Po is the zero potential distance r. If all the mass considered is 
confined to a certain region it may be thought to be concentrated at a point in that region. The mass 

25 may then be considered a point source, and Jv f dV = K = const is independent of the limits of 

integration. In such a case, ris commonly chosen to be infinite, such as in gravity problems since the 
gravitational potential reaches zero at infinity. However, in continuum physics that option is not 
possible since mass and potential are proportional, and k oc V. The potential does not reach a limit as 
V is varied. Since mass is evenly distributed over the region enclosed by the thermodynamic system 

30 and in the immediate surrounding, it represents a distributed source (Kellogg 1929, p. 156). The 
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potentials of distributed source problems are commonly logarithmic. In such a case the zero potential 
distance may be finite, it is then by convention often assigned unit magnitude. 

In thermodynamics, the properties of substances are given in a standard state (P*, V*. r*) 
5 which may then also serve as a zero potential state, and a change of the pressure to ± oo requires 
infinite work. The distance term r is here interpreted to be the radius Tq of the thermodynamic system 
in the standard state. It is the one-dimensional equivalent of the volume of the mass n to which the 
equation of state PV = nRT is scaled. For solids which have a non-zero density in a vacuum, the zero 
potential state is then defined by P* = Pint. A deviation from Tq indicates a change of state. (Cauchy's 
10 continuity approach amounts to the attempt to perform a limit operation with respect to Tq at constant 
conditions which is against the nature of the potential; Kellogg 1929, p.63). In Part I of this series of 
essays, Tq is set to be invariant with respect to direction, i.e. the shape of the thermodynamic system 
is assumed to be a sphere because it minimizes the surface-volume ratio. The material is therefore, 
defined to have isotropic properties. In Part II, additional constraints require another shape. 

15 

Consider the divergence theorem for an isotropic state, such that only normal (radius-parallel) 
forces need to be taken into account, 

jfndA=jV'rdV = fc, (eqn.02) 
where n is a unit vector normal to A. Consider forces as an energy flux, f = eidUldx^ The 

20 surface A is closed to envelop a volume element \dV\ the forces are exerted by the mass in the 
element against the surrounding, or vice versa, ^t similar extemal conditions, a small and a large 
. quantity of mass will do similar amounts of work upon their surrounding, relative to their mass. On this 
requirement of scale independence the entire theory of thermodynamics is based. Thus the : 
thermodynamic system is a source of forces. The divergence of the force field exerted by the system 

25 (and thus also the divergence of the force field acting upon the system) is proportional to mass, hence 
the RHS is a linear function of V; thus div f is a constant. From this follows: at given external 
conditions, div f is insensitive to scale. If the scale of consideration is varied which is measured in V, 
and since the relation of \/oc to >^ oc is not linear, the relation of A to |f| cannot be constant; |f| 
must necessarily be a linear function of scale, |f | oc r. div f is the trace of the tensor F defined below. 

30 Both from eqn.02 and eqn.06 it follows that 
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\-j = div f = const (eqn.03) 

In 

at constant external conditions (Koenemann 2001a), where r is the radius of the system, and 
r is a nneasure of the scale considered. 

These conclusions are not changed by further developments of the divergence concept that 
can only be elaborated below, after explaining some context. Eqn.03 is In agreement with the 
properties of a thermodynamic continuum, with the spatial properties of a thermodynamic system that 
interacts with its surrounding, and with the fact that the system represents a potential of distributed 
matter (Kellogg 1 929. p. 1 56). 

A limit operation with respect to V would not change the relation in eqn.03, but it would vanish 
identically if r reaches 0. This is in accordance with potential theory (Kellogg 1929, p.147), and a ; 
refutation of the Cauchy lemma -f^^ = for which to be valid f must reach a finite value as \/ and r 
vanish (Truesdell & Toupin 1960; Gurtin 1982; Truesdell 1992). 

A thermodynamic system is defined by a chosen amount of mass. Its location in space is 
given by its center of mass Q in terms of external coordinates X/. The surface points P of the system 
are given in the internal coordinates X/ whose origin coincides with the center of mass of the chosen 
system. The physical conditions at the points P In X/ are functions of Q{Xi) in the sense of the 
definition of a vector field. T(Q)v = x. where T is a tensor as a function of location Q(X), v is the 
location vector of a nearby point P relative to Q, and x is a vector located at P as a function of T and 

V. 

The correlation between a surface point P and a particular direction vector in x, is unique. 
Another direction vector passing through P can only be part of a different system with its own origin 
and coordinate set, and different boundary conditions may apply to it if the gradients in X, are non- 
zero. It is therefore not of interest. 

Derivation of the external and material force fields 
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The definition of a tensor is the derivative of a vector field with respect to the coordinates, or 
the second derivative of a scalar field with respect to the coordinates (Malvern 1969. p.57). Energy or 
work are such scalers. Let the external energy be L/ext and the material internal energy be E. so 

^ = /. ^ = m, (eqn.04) 

dx. dx^ 

- =eiF e7 = — ^- = e/M ey (eqn.05) 



dx^dxj dxj dXfdXj dxj 

e/ f = Fp = f ei f ^^%ix, = Mp = m (eqn.06) 

where r is the radius, or the position vector of a point on the surface of the system, f is the 
external force field; m is the material force field, or just the material field; and F and M are the tensors 
10 controlling the properties of the vector fields. F represents the external boundary conditions, and M 
represents the properties of the material which may be underistood as a set of internal boundary - -^^ 
conditions. Forces are understood as an energy flux that changes the internal state of the system, not 
as an external acceleration. 



15 Since externally unbalanced forces cannot cause a deformation, they can be ignored. 

External equilibrium (Newton's equilibrium) is therefore a precondition. Also Ignored are body forces . 
as they interact with the inertial mass, but not with the thermodynamic mass. The internal 
(thermodynamic) equilibrium condition is then given by , 

f + m = 0 . (eqn.07) 
at any point P on the surface of the system, or, If f and m are understood as functions of directions 0, 



dQ = ~jm dQ (eqn.08) 



25 as a sum around the system in 2D. 

The external torque is balanced by definition since the interface between system and 
surrounding is bonded; the disequilibrium case cannot occur in an elastic medium as long as no 
bonds are broken. The condition 



Jf X r i/e = 0 (eqn.09) 
is an equilibrium condition with surprising freedom because f and rmay vary In a reciprocal 

way without changing the result. If dr/dQ = 0, eqn.09 is a statement of orthogonality, describing the 
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properties of the external boundary conditions. The torque of the external forces may be balanced 
with additional help from surface bonding forces nis- The complete condition for the balance of torque 
is 

J(f X r - m^xr)dQ = 0 (eqn.10) 

5 

where tn^ ± r. their magnitude is invariant, 

-!-^ = 0. (eqn.11) 

ae 

10 and ms X r has the same sign at all surface points P. Surface bonding forces are neither 

external forces nor material forces; they are constraint forces which do not do work on either system 
or surrounding, but they make the interaction of system and surrounding possible. Their existence is 
concluded from .the precondition that equilibrium must exist; they balance the torque of f if necessary. 
The subject will .be illustrated later. - 

15 : ^ :^ 

Eqns.07 and 08 are important for the method because they clearly distinguish material force 
and external force, i.e. a material force is recognized as a physical entity of its own. It is exerted by 
the system upon the surrounding due to a change of state in the system resulting from the action of 
external forces. The system thus represents a potential. Equilibrium between system and surrounding 

20 implies that for Isotropic conditions, 

div f + div m = 0; (eqn.12) 
eqn.12 therefore consists of two Poisson equations. 
25 Stress in the sense of this method is a force vector field ftotai resulting from the interaction of 

the two independent force fields f and m, the exterior and the material force field, and their respective 
boundary conditions In the state of equilibrium, plus the surface bonding forces if necessary, 

ftotai = ™> ™s) • {eqn.13) 
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Definition of a materiai-independent equation of state for soiids 
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The ideal gas law, PV = nRT, disregards the atomic structure of matter and considers bulk 
behavior only. Boyle's law can be understood as a material law for ideal gases since it predicts a 
particular behavior of the gas upon a change of external conditions. It relates the internal energy of a 



14 



system, its mass and volume to one another and thus fully describes the energetic state in which the 
system is. In principle, such a description of state must also exist for solids. It must therefore be 
possible to define an ideal solid. 



5 The most general form of the equation of state for a solid is here chosen to be 

P*K = z (eqn.14) 

where 

k= r! ■ (eqn.15) 
The formula predicts that all solids have the same compressibility dVldP if the molar volume 
10 is taken into account, z = /(P) is a number characteristic for a particular state, the function Is not 

known. The quantum-mechanic problem as to how the volume of a solid comes about, is unsolved. It 
is therefore not possible to predict z, but it can be modelled through the Birch-Murnaghan equation. 
The latter is phenomenological. but it is successfully applied in studies of material behavior under high 
pressure, e.g., the prediction of the elastic properties of the Earth's core (Jeanioz & Knittle 1986). 

15 

If the alkali metals are compressed up to'100 kb it is found that they all follow the same 
pattern to a first approximation if the externally applied load is normalized with respect to the internal 
pressure of the solid (Fig.1 ). Thus it seems justified to use the internal pressure as a standard. for the 
behavior of a particular material. Since eqn.14 is Isotropic, k is independent of directions and only 
20 important for the modelling of real materials. For the purposes of this paper k is assumed to be unity 
so that z = const. In that sense, the thermodynamic properties of the ideal solid are those of an ideal 
gas. and Boyle's law applies. 



25 Transformation of the equation of state from scalar form into vector form 



The most basic material law in thermodynamics is the Boyle-Mariotte law, 

PV=^ const (eqn.16) 
The LHS is only the most common interpretation of the product. It may also be interpreted as 

30 a vector form of Boyle's law. 

r * f = const (eqn . 1 7) 
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where r is the position vector of a point P on the surface of the system relative to its center of 
mass. The star product r«f is interpreted as 



must be invariant with respect to direction. r*f gives the work done by an external force f, both 
normal and shear component, at the point P with position vector r on the surface of a thermodynamic 
system. Eqn.18 is known to be an identity. If r is a unit vector, r*f = |f|; however, the ratio of |r| to |f| 

may now be a function of direction while eqn.17 is still observed if 1 f I and | r | maintain a reciprocal 
relation. 

The force field ftotai that develops in the transition from the unloaded to the loaded state can 
be decomposed into an isotropic component, the operative force field fop, and the deviatoric 
component fdev; fop represents the change of state, and fdev gives direction properties. Both f r and 
I f X r I integrated over the surface of the system yield zero. Through fop a reference state is defined 
which serves to assign signs to the directions of the deviatoric field (Fig.03). 

In contrast to gases and fluids, solids can support shear forces, and the energetic state of a 
system is changed by work done by shear forces and normal forces alike. A normal force will cause 
an expansion or a contraction of the system, depending on its sign. The sign of shear forces indicates 
which way the body would spin if they were unbalanced, but it gives no hints regarding volume 
effects. However, the volume effect of shear forces can be only dilational. A deviatoric force field fdev 
can be partitioned into a normal component field fn and a shear component field fs. Assume a 
spherical volume being subjected to the shear components fs of a force field with orthogonal eigen- 
directions, with the origin of the coordinates at the center of mass. Along the surface, all fs will act on 
the surface points P with position vectors r; the points P will be displaced parallel to the direction of fs 
to the deformed position P' with position vector r'. Since r and v = P P' are mutually perpendicular, 
|r'| will be larger than |r|. Therefore the work done by a shear force has a volume effect, and it is 
always dilational. The.dilational component accumulates from zero at the contracting eigendirection 
towards the extending eigendirection where it reaches its maximum. This effect holds both for shear 
forces exerted by the surrounding on the system, and for shear forces exerted by the system on the 




(eqn.18) 
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surrounding. The combined effect of the work done by shear forces results in two additional force field 
components m$(jnt) and fs(ext) parallel to the extending eigendirection, such that due to the existence of 
shear forces a dilational (i.e. normal) effect is produced. Thus 

f total = fop + fdev + m^hear + fshear (eqn.1 9) 

where f = fop + fdev. Fop = cl, and det Fdev = ± 1 . An example is given below. 
Work is then calculated in analogy to PdV-work, 

jfdr = cj— = clnr ' (eqn.20) 



The relation of material distance, i.e. the radius ro, to a collnear force is derived through 
differentiation of eqn.1 8, 

15 . r^df + f^dr = 0 (eqn.21 ) 

If eqn.21 is divided by ro, integrated, and divided by /b. it follows that 

In — (eqn.22) 

20 

in complete analogy to In (V^/Vq) = - AP/Po in isotropic thermodynamics. Eqn.22 provides the 
cause-effect relation, by which the displacement field Ar = s Is generated from the applied force field f 
which therefore have identical geometrical properties. 

. 25 c A solidiin a vacuum has its ideal volume Vq with unit radius Tq. The internal pressure Pint of a 

solid is in the order of several kbar. In directional terms, the analogue to Pint is the internal force mo 
which is a non-zero dormant force as it is internally balanced In the unloaded state. For modelling 
purposes it is set to unity in the following text. Above it was explained that the work done by shear 
forces on a volume is a dilation. Thus, be it a normal component fn or a shear component fs. their 

30 combined effect, colinear with X/ needs to be considered. The magnitude of the externally effective 
force must be scaled to that of the material force mo, so eqn.6 in complete expansion is 

-Am = Af.,. = 1^1^^ (eqn.23) 
l™o| |fcx.| 
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where the first RHS term gives the magnitude of fexi in multiples of the internal pressure of the 
material, and the second term is a unit vector with the orientation of text- This normalization is always 
implied; subsequently unit magnitude is used, and Af is simply referred to as f. 

The condition of equilibrium for the torque is dependent on the properties of the force field 
and the shape of the system. The variability of the shape is constrained by the properties of the 
external force field, the material properties, and the condition that system and surrounding are 
bonded. If the effect of a progressive deformation is demonstrated in 2D, commonly a circle is trans- 
formed by a displacement field into an ellipse. The unit circle of the undeformed state as a geometric 
device is certainly adequate for this purpose; but the shape of the thermodynamic system (or volume 
element) must fulfil the equilibrium conditions. Specifically, the thermodynamic system may have an 
elliptical shape in one particular set of conditions, but It may be a circle in another. The two concepts - 
a circular pattem of points in the undeformed state vs. the mechanlcially active shape of the systenfi - 
must not be mixed up. ^ 

Strain is a function of the displacement field. If the principal axes of the strain ellipsoid do not 
rotate during progressive deformation they must coincide with the characteristic directions of the 
displacement field (and thus the force field). It can therefore be concluded that the eigendirections of 
the force field are mutually perpendicular, or that the force field, and hence the displacement field, are 
orthogonal. From a mathematical point of view this is a rather special case, however, and by no 
means a precondition. The strain tensor is thus not a helpful term to understand the physics of 
defonmation. It is therefore necessary (a) to establish the shape of the system from the equilibrium 
conditions, and (b) to find the eigendirections of the force field. They are defined as the directions 
along which the force field f has no shear force components. 

Shape of the volume element 

Eqn.18 gives constraints for the shape of the system. The simplest shape is the sphere - it is 
isotropic, minimizes the ratio, and for any point source at the center of mass its surface is an 
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equipotential surface. For isotropic external conditions and an isotropic material a spherical shape for 
the system is the most natural choice as there are no shear forces. For other deformation types 
eqn.18 Implies that radius and normal force along the contracting eigendirectlon c and those along the 
extending eigendirectlon e are identical in magnitude. The curl of the pure shear deviatoric field 
(Fig.04) 

" ¥2 



fix. 



0 

dx 



1 J 



(eqn.24) 



is zero; therefore the eigendirections of the resulting force field coincide with those of Fdev 
The conditions tr T = 0 and det T = ± 1 are both conservation conditions und must hold 
simultaneously, hence I Fn I = I F22 1=1. Since F<jev is orthogonal, there are no constraints from the 
equilibrium conditions on the shape of the volume element; it is therefore still spherical in shape. 



The magnitudes of normal and shear deviatoric forces are shown in Fig.05. Under the given 
boundary conditions, all deviatoric forces on the surface have the components f = [cos9, -sin9] where 
15 e is the angle between the orientation of r and the Xi-coordinate. The external equilibrium condition for 
the torque for the entire body is then . 

xr Je = 2j^"cose sine Je = 0. (eqn.25) 

i.e. it is balanced without the help of the surface bonding forces nis. • 

20 

Kinematics of pure shear 

mo is the material force which is in equilibrium with a vacuum. For a body with unit radius in 
25 an orthorhombic deviatoric field, the components of the position vector r of surface point P are [r^ rJi = 
[cosO sinG]; the operative force vector is a normal force with components fop proportional to [-cos9 
-sinG]; the deviatoric force fdev is proportional to [cosG -sinG]; and the radius-normal unit vector t = [sinG 
-cosG]. 
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External forces acting on point P are the isotropic operative force fop and the deviatoric 



component fdev fop consists of normal components only, fdev can be split into a normal component 



^n(dev) and a tangential component fs(dev)- The sum of all normal components is 



2f „ = f,p + (r . f )r = 1 + (cos^ 9 - sin^ 9) 



(eqn.26) 
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The respective inward-directed displacement vector Sn parallel to r is found through eqn.22. 



The cross product 



»• X fdev I = |f s(dev) I = 2 cos9 sin9 



(eqn.27) 
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is equivalent to f • t . A shear force has the effect that P is displaced away from the 
contracting eigendirection, and towards the extending eigendirection (Fig.06). The components of 
fs(dev) at P evoke elastic material forces of equal magnitude and opposite sign within the system such 
that the torque at P is balanced. The effect of fs is that P is displaced to P' (Fig.06); however, since fs 
15 is perpendicular to r, the distance from the origin to P' will always be larger than that from the origin 
to P (Fig.06). That is, a shear force will always cause a local extension proportional to its magnitude. 
That effect is additive from point to point. The observed dilation at P* is therefore the sum of all dila- 
tions from the contracting eigendirection to P, 



where a is the angular distance measured from the contracting eigendirection c at X2 to P. For 
a = "/2, Jf t d9 = 1 . At the lower limit of integration there is no dilational effect, naturally; the integral 
reaches its highest value at the extending eigendirection e although no shear forces are observed at 
that point. 



In a continuum.'the result of eqn.28 counts twice. Assume a body loaded hydrostatically so 
that the operative force field (the change of the energetic state of the system) has magnitude -1 , and 

30 the deviatoric field has magnitude zero. Along Xi and X2 the radius is contracted. Then the boundary 
conditions are relaxed along only; consider the surrounding as immobile. In the attempt to reduce 
the stored energy, the material expands by itself parallel to Xi whereas nothing changes along X2. The 
operative field magnitude is reduced to -14. and a deviatoric field of magnitude ± 14 develops. Since 
the expansion is not isotropic, the body will exert shear forces at the surrounding; in the first quadrant 

35 (to the right of the contracting eigendirection and above the extending eigendirection) it will result in a 




(eqn.28) 
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Figure 5 here 
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sinistral couple at the system-surrounding interface, and in a dextral couple in the second quadrant 
(counting clockwise). Therefore, work is done by the system on the surrounding, thereby reducing the 
energetic state of the system from -1 to until the system has reached its energetic minimum state. 
Note, however, that the sign of the couples is opposite to what should be expected from a proper 
5 continuum pure shear case. But if the surrounding deforms with the system. It will first expand with the 
system in x^, and then exert additional shear forces at the system, causing a dextral couple in the first 
quadrant and a sinistral couple in the second, resulting in a further expansion parallel to X2. This shear 
work is work done by the surrounding, and will cause an Increase of the energetic state from -V4 to -1 . 
Therefore one part of the expansion is supplied by the shear effect of the material, one part by the 
10 shear effect of the external force field, and both are proportional to the magnitude of the devlatoric 

field. The result is that a dilational force of twice the magnitude of the integrated external shear forces 
acts upon the surface point on the extending eigendirection which will therefore be shifted away from 
the center of the. system. 

15 

Volume effect ofdeviatoric loading 

A numerical example, using eqn.22 (inward-directed work is positive): Tq = 1, /b = 1, and 
volume Vo/n = 1 . After loading to operative field magnitude 1/5, the magnitudes along X2 are Af = fop + 
20 fn(dev) = 2/5, r2 = 0.670, and along Xi they are fop + /n(dev) = 1/5 - 1/5 = 0, ri = 1 .000; but /s(int) + 4(ext) = - 
2/5, so Ti = 1 .492. The volume of the resulting ellipse is V/n = r^r2 = 1 , i.e. area is preserved if 
compared with the unloaded state (Fig.07). 

25 Model calculations for simple shear 

The eigendirections 

The position O of the body center in space is given in coordinates X/. A second coordinate set 
30 X/ is fixed to the body, its surface points are given relative to O in the latter system. 



21 



For a dextral shear the deviatoric external force field therefore has the form Fextr = fexi where 



F = 

* ext 



0 0 



(eqn.29) 
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(Fig.OSa). This field interacts with a unit body (the thermodynamic system) possessing 
isotropic material properties (Fig.OSb). Thus, in a continuum an infinite number offeree vectors will act 
at an infinite number of points of action along the surface of the system. Therefore it is necessary to 
find the average point of action of the average force vector fav on the body. From the statement of the 
problem (Flg.08)i foe sin 9 where 0 is the angle enclosed by r and the reference line Xi. The average 
for one quadrant is found through 



^sinG^/G 2 

71 



n 
2 



(eqn.30) 



The average force is a vector of the form fgv = [2/7t, 0], and the position vector of its point of 
action encloses with x^ the angle k = ± 39.54'' (Fig.09). If the body is fixed in space and subjected 
to a force field as in eqn.30, the normal force component has the form 

f „ = (f • n)ii = cos G sin G , (eqn.31 ) 

and the shear force component is given by 



At the latter is therefore 



=(f .t)t = sin'G. 



f2y 

fs(K)= - =0.405, 



(eqn.32) 



(eqn.33) 



If the body is allowed to react to the applied force field, fs(K) is accommodated: since the 
coherence is maintained, external disequilibrium is impossible. Hence the system is not able to rotate 
freely, but surface bonding forces ms will be activated parallel to its surface. Thus |fs{K}| is subtracted 
from all |fs| at all points; the sign of the effective fg will be reversed in some areas. The effective shear 
force magnitude at any point P is therefore 

|fs(eff)|e=fe-te-|fs(K)|. (eqn.34) 
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and the total effective force at P{B) is 

fefr=fs(c«)+fn- (eqn.35) 

5 

Initially there are two points with no radius-normal force components at 9 = 0 and Q = n 
(Fig.OSb). Due to the subtraction of fs(K). these points move away from Xi by 0 = ±-k to either side of 
the coordinate axis. Since only normal components act along these directions, they are interpreted as 
eigendirections (contracting: c, extending: e) which are not mutually perpendicular. They enclose the 

10 angles ti - 2k = 100.92** and 2k = 79.08°. However, by using the unit vector n and t ± n in eqns.31 
and 32. a spherical shape of the system is implied because there is no difference yet between n and 
r. It is still not possible to balance the rotational momentunri for a body with a spherical shape 
as I fs(eff)^r | dO^O. Therefore this assumption regarding the shape may be faulty. The angular 
relation of the eigendirections suggests that the effective force field fetr for simple shear and the shape 

15 of the system both have elliptical properties. It is therefore necessary to calculate the shape of the' 
ellipse which is in equilibrium with feff with eigendirections as indicated above, and which fulfils the 
requirement of eqn. 18. 

20 Force field and shape of volume element 

In the pure shear example (PS), the field property tensor (eqn.24) has the properties 
tr Fps = 0 and det Fps = ± 1 . Eqn.29 for simple shear (SS) differs in the latter point. However, the term 
in eqn.33 is not part of eqn.29. The field matrix for simple shear is therefore assumed to be 
25 characterized by the condition det Fss = ± 1 after the subtraction, indicating that all normal fluxes that 
are part of the deviatoric field and which enter the system, will leave it as well. Thus, the force field 
under the boundary conditions for elastic simple shear and the shape properties of the system are 
here modelled by assuming an orthogonal field of the form 




(eqn.36) 
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which is in equilibrium with a system of spherical shape. The maximum shear directions 
should coincide with the points on the coordinates, though, so the force field is reoriented by a 
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Sinistral rotation of 45°. Both force field and radius field are then transformed by a transformation 
matrix T to acquire elliptical properties, where p are the position vectors of the surface points in the 
untransformed state: 



rcotK 0 

[ 0 tanK 
[tan 



(eqn.37) 



COtK 

K 0 

Tp = = [r, , cos 9 7^22 sin 0] , 
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and 



{eqn.38) 
(eqn.39) 

FdevP = fdcv =[7;,sine cos 0]. (eqn.40) 
The eigendirections of the elliptic radius field ten are thus mutually perpendicular, and identical 
to the principal axes of the resulting ellipse with area A/n = TnTaa = 1 whereas the eigendirections - 
extending: v©, contracting: Vc - of the force field are non-orthogonal as desired. 



' 1. 

± tan^J 



(eqn.41) 



Kinematics of simpfe shear 

Conveniently, all equations regarding equilibrium conditions etc. are just transformed forms of 
the pure shear example, too. The cross product simplifies to 

^Ifdev X nil 1^0 = T,,T22 ^(cos' 9 - sin' 0^9 = 0. (eqn.42) 
i.e. the elliptic properties of radius field and force field cancel. The linear momentum is 

jf^^^ • r^„^/9 =(7;^ + 722 j^cos 9 sin QdQ = 0 . (eqn.43) 

Since T^] + > 2 , the maximum magnitude of the dot product along the eigendirections is 

larger than unity. This Is believed to be an artefact of the elliptic shape as the radii along the 
eigendirections do not have unit length due to the transformation by T. If the force is dotted with the 



inverse position vector Ten"^ = tV. 



^|f„(dev)|^6 =j^cn -fdcv^Q =^2cos0sin0i/0 = 0 



(eqn.44) 



the equilibrium condition Is still maintained. The elliptic properties of force field and body 
shape cancel, resulting in a consideration of normal deviatoric force f(n)dev per unit radius, the 
maximum magnitude of which is 1 . Eqn.44 gives the normalized relation of |fn(dev)| to |reii|; that 
normalization of fn with respect to r is necessary (cf. eqn.23). fop (cf. eqns.19 to 23) is therefore also a 
unit vector. Eqns.42 and 44 together observe eqn.18. Some vector magnitudes are shown in Fig.1 1 . 
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If t ± n, t is not a tangent vector to the ellipse, but is defined as a radius-nomrial unit vector. 
Irrespective of the orientation of the elliptic surface. 

5 The magnitude of the shear force component is given by 

|fs(dev)| = *^dev • t = ^,1 sin' e - cos' 9 (eqn.45) 
Integration for separate sectors (e.g., from the contracting eigendirection to either side 
towards the extending eigendirection, or from 0 to k and from k to nl2) shows that there is an • 
10 imbalance between dextral and sinistral shear forces acting on the system. However, eqn.45 is not an 
equilibrium condition, and the shear forces by themselves need not balance. Because forces are 
balanced for the ellipse (eqn.43), an imbalance in fs(dev) cannot imply a spin. The imbalance is 
therefore interpreted to indicate a permanent reorientation of the force field through a transformation 
of the internal coordinate axes X/ with respect to the external reference frame X/ by an angle X such 
15 that 

V|fs(dcxtra,)|-|fs(sinistrai)| = tBH^l; X = 28,83^ (eqn.46) 

where the vector magnitude terms under the root stand for the integrated sums over adjacent 
sectors. Since the eigendirections are real, the transformation angle X is stable. The root is taken 
20 because all angular terms in eqn.45 are squared. 

It is useful to interpret the LHS of eqn.46 as a tangent term because there is no reason why it 
should be restricted to values ^ 1 . As the LHS oo, A. -> QO"*, which would be the condition for 
external rotation. At the same instant, the characteristic equation should degenerate, whence the 
25 eigendirections become imaginary. 

The invention contains the notion that a shear force has the effect that P is displaced towards 
P*. away from c, and towards the e. A shear force will always cause a local extension proportional to 
its magnitude. That effect is additive from point to point. During a simple shear elastic deformation, the 
30 same principle holds, except that the eigendirections are no longer mutually perpendicular. However, 
the result is similar to that of eqn.28 (which was 1 over each quadrant), except that the surface over 
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which is the integration is perfomned, is larger than that of the circular volume element. The result per 
quadrant is therefore larger than for a circular surface by the factor TuTzJl = 1 .018. As explained in 
the pure shear example (following eqn.28), this result counts twice, = 2.036. 



5 Ws is the dilational effect of the shear forces along the radius parallel to e. The shortening 

factor for the radius Tc along c is given by eqn.22. If Vc and are the normalized eigenvectors, they 
can be understood as the radii of a unit circle subjected to the monoclinic field Fr = f , and the 
displacement field can be calculated. A numerical example, using eqn.22: Tq = 1, fo = 1. and volume 
Vq/tt = 1. After loading to operative field magnitude 1/5, the magnitudes along c are = 

10 2/5, Tc = I vj = 0.670. Along e they are + /n(dev) = 1/5 - 1/5 = 0, Te = 1 .000; but Ws = -2.036/5 = -0.407, 

^ so Te = I vj = 1 .503. If S is the displacement field property tensor^ 

' r^r^ =|vJ|vJ = detS =1,0075 (eqn.47) 

indicating that in this example the volume expands by 0.75%. 

15 

Geometric properties of simple shear deformatior) 

The coordinates are rotated with respect to X. by X in the sense of shear. In physical spade 

20 X/, the following picture emerges (Fig, 12): the extending eigendirection e is at 10.71"*, the contracting 
eigendirection c is at 1 1 1 .63". The bisectors of the angles enclosed by c and e are parallel to the Xy 
coordinates and represent maximum shear directions. The orientation vector bisecting the large 
sector at 61 ,18** indicates a shear plane perpendicular to that orientation, the R-plane is therefore 
found at -28.83*'; the orientation vector bisecting the small sector at -28.83** indicates the P-shear 

25 plane at 61.18** (Fig,12). 

The topic of this invention is the determination of the force field that causes the deformation, 
and the displacement field. But irrespective of the deformation mode - elastic or plastic - the 
displacement field must reflect the properties of the force field even if the resulting features may differ 
30 considerably in their nature. Therefore it is justified to correlate the properties of the calculated force 
field with fabric elements observed in plastically deformed rocks. X, is the shear zone boundary 
(Fig. 12a). Of the eigendirections, e is identified as the S-plane in S-C-fabrics (Berthe et al. 1979) and 
the plane parallel to which the main anisotropy of crystals is oriented. It is thus a stretch-only, no- 
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rotation-no-shear direction. The predicted angle of ca. 1 1 ** compares favorably with the obliquity of 
fabric diagrams from monomineratic shear zones formed of minerals with only one major shear plane 
(ice, mica, olivine). 

5 c is not well developed in high-temperature tectonites. It may be recognised through a lack of 

pressure shadows, or minimum mica alignment along the surface of feldspar porphyroclasts. . 
However, if mylonites are exhumed they commonly develop joints that cut the layering at ca.70-80°, 
consistently inclined against the direction of shear. They appear to be controlled by the elastic energy 
stored during plastic deformation which is released when the confining pressure is no longer able to 
10 hold the rock together. This invention here predicts that dilational cracks should open parallel toe 

The shear plane at 9 = -28.83* is the Riedel plane R (Riedel 1929) or the C-plane in S-C- 
fabrics (Berths et al. 1979). The P-plane is usually suppressed in natural plastic defomnation, but has 
been observed in shear box experiments (e.g., Merzer & Freund 1975), and is occasionally found as a 

15 minor shear direction (Fig. 12b). It is better developed if the rock behaviour was near the brittle-plastic 
transition zone. The synthetic R-plane (here dextral) and the antithetic P-plane (here sinistral) are not 
mechanically equivalent. The R-plane is expected to be the first to yield at the onset of plastic flow. 
Fig.13 shows the flow field. The sense of shear on the R-plane is synthetic to the bulk sense of shear, 
but the R-plane.also simultaneously rotates antithetically towards e; in the process it is progressively 

20 stretched until it decays. R- and P-planes are stretching faults in the sense of Means (1990). 



Energetics of elastic and plastic, pure and simple shear defonnation 

25 

Elastic deformation 

The state in which the system is, is controlled by the operative force field which serves as a 
measure of the field strength. It can be used to compare deformation types with one another. In the 
30 elastic mode, the operative and the deviatoric field are proportional to one another during loading; 
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thus they both require work which represents the elastic potential. Since the sum of fn(dev) from 0 to h 
is zero, the work done by normal forces is only the work done by fop, 

f/op = -271 (eqn-48) 

5 for To = fo = Afop = 1 (cf.eqn.22). For an isotropic state, this would be the only term to be 

considered in 2D; in absolute numbers it tums out the be the maximum work for all states. This work 
will cause an isotropic contraction. Work done by shear forces will cause a directed dilation. 
Integrating eqn.28 from 0 to 7i/2 for one quadrant, times 4 to account for all quadrants, yields 8.000; 
total work in elastic pure shear is therefore wps(ei) = 8 - 27i = 1 .717. 

10 

In elastic simple shear, the magnitude of fop is the same as above. The shear work term 
(eqn.47) must be evaluated for one sector from c to e, and multiplied by 4 to account for all sectors. 
Total work in elastic simple shear is thus Wss(ei) = -27t + 8.144 = 1 .861 . At similar operative field 
strength, an elastic simple shear thus requires 8.4 % more energy than pure shear. 

15 

Piastic deformation 

The mechanism by which plastic deformation is achieved is not the subject of this invention; 
20 however, it is here assumed that the material remains continuous, and that a homogeneous state of 
the force field f is maintained. The transition from the elastic to the plastic deformation mode (the yield 
point) is marked by a change from reversible to irreversible behavior. Thus the operative field reaches 
a stable value requiring no more work. 

25 All further work, the plastic deformation energy, is then work done by shear forces only. Since 

the dilation caused by them is also an elastic feature the buildup of which stops at the yield point, the 
energy is not stored in the volume of the system, but dissipated along its surface. Since the system 
therefore behaves passively, the double shear work concept (see discussion following eqn.28) no 
longer applies, thus the required shear energy is only half the amount for elastic deformation. 

30 Dissipated energy is energy lost and need not be balanced; however, since no free whole-body 
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rotations are possible In a field with real eigendirections - which are still required by the continuing 
existence of the elastically loaded state - the rotational momentum must still be balanced. 

For pure shear, the sinistral and dextral shear integrated balance, thus no external rotation of 
the system is possible. In terms of work done, sinistral and dextral shear work add, thus the work 
done in plastic pure shear is half of the amount done in elastic loading, i.e. wps(pi) = 4. 

For simple shear the affairs are a little more complex. Eqn.45 integrated from the contracting 
to the extending eigendirection for both sectors yields differing magnitudes: 1 .339 dextral shear work 
in a large sector, and 0.733 sinistral shear work in a small sector. The Imbalance for the entire system 
is 1 .202, the total shear work done on the entire system is 4.144. Since the large sectors dominiate, 
the imbalance is resolved into plastic dextral bulk shear. However, the resolution of the dextral 
imbalance leads to an imbalance in the angular momentum of the entire system. Since the shape of 
the system is elliptic and the long axis of the system is in the small sector with sinistral shear, the - 
angular momentum of the system is proportional in magnitude, but opposite in sign to that of the 
resolved shear, i.e. sinistral. The dextral plastic shear work costs energy which is dissipated, i.e. it is 
an internal dextral rotation. The rotation of the system due to the angular momentum is an extemal 
sinistral rotation which is free. Since internal and external rotation cancel, the eigendirections are still 
real. Thus the energy dissipated due to a plastic dextral shear is minimized by an external sinistral 
rotation, and the energy thus saved is proportional to the imbalance; thus Wss-pi = 4.144 - 1 .202 
= 2.933. 

At similar operative field strength, a plastic simple shear thus requires 26.7 % less energy 
than pure shear. Energetic differences of such a magnitude have been observed in experiments; they 
occurred systematically as a function of the displacement field, and independent of the experimental 
substance (quartz, salt, gypsum; Franssen & Spiers 1990; Franssen, pers. comm.). 

The procedure presented here is completely compatible with the principles of potential theory 
and thermodynamics. The invention offers a straightforward, immediate correlation of force field to 
displacement field which are both vector fields controlled by tensors; from a mathematical point of 
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view this is more systematic than the procedure used so far. The new method also permits predictions 
which are not offered by the known procedure because the latter does not permit e.g. the deduction of 
a particular displacement field from a given stress state. 

5 A viscous deformation step for a gas or fluid can be decomposed into a time-independent, 

reversible, elastic loading step and a time-dependent, irreversible, diffusion-controlled relaxation step. 
The new method is also suited for the modelling of viscous fluids as long as time spans are 
considered which are short in relation to the relaxation time. 

10 . 

An example for the modelling of an elastic deformation of a sample with given shape under 
given external boundary conditions 

The approach applied here is based on potential theory. The thermodynamic continuum is not 
15 a continuum of points, but a continuum of systems. Any point Q in Euclidean space within a region of 
distributed mass can be thought to be the center of mass of a thermodynamic system; its physical 
properties - density, chemical composition, temperature, material properties, orientation of anisotropy 
- plus, those of the external boundary conditions are then functions of Q. Scale independence ensures 
that the extensive properties are scaled per unit mass. Thus two infinitesimally close points Qi and Q2 
20 represent two different thermodynamic systems V\ and V2 of some finite unit size; they largely 

overlap, but their physical properties or boundary conditions may be subject to gradients in Euclidean 
space. The systems may be understood to have unit mass, unit volume, or unit radius which are all 
finite. In contrast to the known method, the invention presented here offers a finite distance which is 
based on its inherent logic - the radius of the thermodynamic system - which can be used as a base 
25 for a wave function. The spatial properties of the physical entities can thus be modelled through an 
application of Fourier series methods as a function of the external boundary conditions. 

The total force field is partitioned into the isotropic operative force field fop and the deviatoric 
force field fdev (cf. Part I, Koenemann 2004b). The latter can again be decomposed into the norma! 
30 component fn(dev) and the shear component fs(dev)- Normal forces will shorten or stretch the radius 
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vector on which they act, depending on sign; shear forces will always have a dilating effect; which 
amounts to an additional stretch component parallel to the extending eigendirection e. Following 
eqn.28 (Koenemann 2000) it was stated that the work done by shear forces counts twice because one 
part of the stretch due to shear dilation in e is done by the system itself, the second part is done by 
5 the surrounding. The component originated by the system is purely a function of the loading 

configuration (loading intensity, external boundary conditions, e.g. pure shear). The shear dilation 
caused by the surrounding, however, surely exists in a perfect continuum, but its effect must reach 
zero at the interface to free space. If all other possible variables are kept constant it is this component 
that must vary as a function of location within a body, i.e. it is sensitive to the body shape and the 
10 specifics of the loading configuration. 

The 2D-model consists of a rectangular body of solid. It is thought to be loaded on top and 
bottom with constant force whereas the lateral faces are free. The applied load is vertical, or parallel 
to y {dfy/dx = 0). The loaded faces (y = 0, y = 1 ) may not change length in x due to friction at the 
15 pistons; at all other y they are free. Thus the load is the same throughout the body, there is no 

gradient. This includes the points on the lateral faces. Consequently the magnitude of the operative 
field fop and the normal force field fn(dev) are monotonous throughout the body. They are not included in 
the model here. ^ 

20 The boundary conditions for the dilatlonal shear 7 vary over the body in a two-fold way: (a) 

the restriction of no change of length in x suppresses even the dilational effect produced by the 
system (system effect, Ti), it is strictest along the. loaded faces y = 0 and y = 1 , and relaxes towards 
the center of the solid, allowing for progressive expansion in x; and (b) the shear effect produced by 
the surrounding (surrounding effect, T2) is fixed to a given minimum, magnitude - on the free faces x = 

25 0 and x = 1 : because there is no bonded surrounding any more; on the loaded faces: because 
expansion in x is prevented by condition (a), but it increases towards the center. The direction of 
relaxation is towards the center for Ti, and towards the free margins for T2. Thus the two boundary 
conditions are opposite in effect and sign (Fig. 14, arrows in fourth ring for conditions 1 and 2). They 
need to be modelled separately. 
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In the static loaded state system and surrounding do work on one another, but we cannot 
measure it; in thermodynamics only the work done in a change of state can be measured. No time- 
dependent processes are observed, system and surrounding are in equilibrium with one another, all 
fluxes are balanced. The quantity under investigation T= + Tz 'is the work done by shear forces, a 
scalar quantity, resulting in a volume change, also a scalar quantity. It is thus possible to approach the 
modelling of the static state through a Laplace equation. 

Ti will be dealt with first. The distribution of a scalar quantity U over a region x(0->d)y(0-»cO 
with a maximum U = M for all points y = 0, and fi = 0 for all points x = 0. x = i?. and y = d is 



mn{d — y)^ 



. rmvc 

sin 



t,=-AMY ^ ^ — (eqn.49) 

- b 

^ (Boas 1983:546. eqn.2.17) where M is a magnitude scaling factor, b is the width of the solid 
along x, d is its length along y, and m = 2 n - 1 (m is always odd). The reverse case is that of ti.= M 
.15 for all points y= of. and /W=;0 for all points. X = 0, x = and yi= 0, . - v 

. ,f mny^ . rmvc 
smh\ — — \sin 

= - 4M y ^ ^ — (eqn.50) 

T:t .1 dmn 
* mnsinh 

b 

The combined effect is Ti = fi + f2. Along x = bl2 forms an upward-warping curve with a 
20 minimum Ti = for y = 0 and y = d, and a maximum {M<0) for y = d/2. Along y = d/2 T^ forms a 

• downward-warping curve with a maximum = 0 for x = 0 and x = b, and a minimum at y = d/2 (which 
is the maximum point on x = 6/2). 

For T2 the boundary conditions are: M = 0 for all points x = 0, x = 6, y = 0. and y = d. and 
25 increasing towards the center. This is done through 



r, =bsm — dsin— (eqn.51) 
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The sum 7= 7, + is interpreted as the potential for dilational cracking (Fig.15-18, upper 
panels). Of special interest is to find out were the gradients of T are the largest. Thus 



f , - mn(d — v) . , mny^ ( mwc 

sinh —+'Sinh — — \cos\ 

,1, b b ) y b 



^ = - AMY^ -J-^ "^-^ (eqn.52) 

^ - bsinh^ 

b 

and 

dT^ Ttx ^„,^ny 



dx 

the magnitude of the gradient is given by 



- n cos — d sin — ; (eqn.53) 



. , mn(d-y) . , mny^ . rmix 
Sinn — + Sinn — ^ \mn sin - 



10 ^ = 4A/t^^— ^ " ^ ^ " ^ (eqn.54) 



and 



2 . to: , , 7t V 
^2rr 71 sin — dsin — 



(eqn.55) 



dx^ b. 

The derivative ^Tldy? (the sum of eqns.54 and 55) of the surface T{b,iS) is interpreted as the 
15 potential for shear cracking (Fig. 15-1 8. lower panels). The derivative d^Tldy^ gives identical results in 
the sense that the morphologies of both derivatives are identical, but their absolute magnitudes are a 
function of the chosen absolute dimensions of b and d, respectively. For realistic results the model 
needs to be rescaled to real material properties. 

20 The method offers all properties of the loaded state as a function of locationlwlthin the region 

of validity bd. In a similar way it should be possible to derive the displacement fields from the work 
equation and to calculate the shape of the deformed body. 



25 Results 



In Figures 15-18, the upper panels show the potential for dilational cracking T, the lower 
panels show the potential for shear cracking ^Tldy?. The width of 10 distance units is the same In all 
figures. Colors indicate relative magnitudes (purple = minimum, red = maximum). 
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For thick, short columns (Fig.15. height 2 units) the center is dynamically dead, it is under 
strong compression. The dilational crack potential is very close to the free surface, suggesting that a 
brittle material loaded to the elastic limit might spallate. The shear crack potential is farther inside the 
5 body, further destabilizing the surface, and low In the center. The failure-potential distribution is 
strongly controlled by the zero-glide boundary condition along the loaded faces. 

If the height of the body is only slighty increased (Fig. 16, height 5 units), the dilational crack 
potential is still close to the surface, but the shear crack potential shifts very quickly into the interior of 
10 the body. The dilational and shear maxima no longer coincide. 

In a body of square shape (Fig.17, height 10 units) the dilational crack potential and the shear 
crack potential again coincide, but this time in the center of the body, suggesting that any material 
yield, brittle or plastic, should start there, 

15 

Turning to a long, slender column (Fig. 18, height 30 units) the dilational crack potential stays 
in the center of the body, but it extends along the long axis, decreasing only close to the loaded faces. 
The potential for shear cracking has left the center and shifted towards the loaded faces, it fonms two 
maxima which peter out towards the center of the body, but they form sharp bifurcating extensions 

20 towards the corners of the body where the loaded and unloaded faces define a discontinuity in the 
boundary conditions. The two extensions and the loaded body surface form a triangular region which 
is also under tow dilational effect; this triangle is largely dead. The calculations permit the 
interpretation that brittle yielding should start through dilational cracking in the center of the body; the 
cracks should propagate along the long axis (which will ultimately change the boundary conditions to 

25 some degree) up to the region where the dilational potential abates; but here the shear potential 
reaches its maximum, so the mode of yield should switch from dilational to shear cracking and 
propagate towards the corners. The material response predicted here is indeed observed in 
experiments on brittle materials such as ceramics (Winter 1992). 
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Conclusions 
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It was shown that the new method for the determination of deformation can be modulated 
such that the physical properties of the loaded state in a discrete body subjected to a specific set of 
5 external boundary conditions can be calculated. In comparison to the known method the requirements- 
on technical equipment and mathematical skill to achieve the results presented here are extremely 
modest. The calculations can be done on standard computers which are widely available. 



The term "computer" is not to be understood in any way as limiting. It can be an arbitrary unit 
10 suited to conduct calculations, e.g. work stations, personal computers, microcomputers, or a switch 
designed to conduct calculations. 

It is a substantial advantage of the method presented here that a grid of pre-defined points is 
not necessary; instead, a solution can readily be found for any point within the region of validity, using 
15 mathematical methods for which solutions for problems in other field theories are known for a long 
time. 

It is of special advantage of this invention to determine deformations without a grid of nodes 
(cf. Fig.01). r 

20 

The invention thus offers a substantial simplification for the determination of deformations. 
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